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Abstract 



We construct new pairs of orthogonal maximal abelian *-subalgebras 
of Me{C), by classifying all self-adjoint complex Hadamard matrices of 
order 6. In particular, we exhibit a non-affine one-parameter family of non- 
equivalent Hadamard matrices of order 6. In the last part of the paper 
we present other previously unknown examples of complex Hadamard 
matrices of higher orders. 

1 Introduction 

Let (Af„(C), Tr) denote the algebra oi n x n complex matrices with the 
usual trace. For X e M„(C), denote by X* the conjugate transpose of X. A 
subalgebra A of M„(C) is called a MASA if it is maximal abelian and closed 
under the * operation. It is easy to see that ^ is a MASA if and only if it 
is unitarily conjugate to the algebra I?„ of n x n diagonal matrices, i.e. A = 
UT>nU* for some U S M„(C) unitary matrix. 

We say that two MASA's Ai,A2 are orthogonal if y^i n .42 = C and the 
vector subspaces Ai QC and © C are orthogonal, with respect to the inner 
product < x,y >= Tr{y*x), x,y G M„(C). This is equivalent to saying that 
the square of inclusions: 



is a commuting square in the sense of [Pol],[Po2] (see also [GHJ]). 

We may assume, up to unitary conjugacy of commuting squares, that Ai = 
Vn and A2 — WDnll* , for some unitary U G M„(C). In this notation, the 
orthogonality of Ai , A2 amounts to U having all entries of the same absolute 
value -j=, hence H = ^JnU is a complex Hadamard matrix. Thus, being given 




Ai C M„(C) 
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a commuting square £ of this form is equivalent to having a complex n x n 
Hadamard matrix. 

Recall that two complex Hadamard matrices are equivalent if there exist 
unitary diagonal matrices Di , D2 and permutation matrices Pi , P2 such that 

H2 = P1D1H1D2P2 

It is easy to see that equivalence of Hadamard matrices corresponds to iso- 
morphism of commuting squares, via the identification described above. Our 
interest in Hadamard matrices, mainly in obtaining one-parameter families of 
non-equivalent Hadamard matrices, comes from the possibility of constructing 
subfactors from the corresponding commuting squares (sec for instance [Jo2]). 
However, it is hard to decide if such subfactors are non-isomorphic, or to com- 
pute their principal graphs. 

Besides their connections to von Neumann algebras ( [H] , [H J] , [.Jo2] , [MW] , [Ni] , 
[Pe] , [Po2] ) , complex Hadamard matrices have numerous other applications such 
as the theory of error correcting codes ([CH]), spectral sets and Fuglede's con- 
jecture ([T]). They play a very important role in quantum information theory, 
in the construction of tclcportation and dense coding schemes ([We]). 

While all Hadamard matrices of orders up to 5 are classified ([H]), it seems 
very hard to describe Hadamard matrices of higher orders, such a classification 
not being known even for n — 6. For n composite, some constructions of para- 
metric families of Hadamard matrices whose entries are linear functions (also 
called affine families, sec [TZ]) arc presented in [Di],[MRS]. There is however 
no general procedure of constructing such families with non-affine entries, or 
for n prime. A catalogue of most known complex Hadamard matrices of small 
order (up to order 16) can be found in [TZ]. 

In this paper we classify, up to equivalence, all Hadamard matrices H of 
order 6 that are self-adjoint, i.e. H = H* , where H* denotes the conjugate 
transpose of H. We thus obtain a new one-parameter non-affine family: 
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where 9 G [—'K,—arcos{ ^^^ )] U [arcos{ ^+^ ),7r] and the variables x,y,z,t 
are given by: 

1 + 2?/ - 



y = exp{iO), z = 



l + 2y + y^- 72^1 + 2y + 2y^ + y^ 
1 + 2y - y2 

1 + 2t/ + - x/2 yi + 2j/ + 2j/3 + 
-l + 2t/ + y2 
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Our family is non-affinc, in particular it is not obtained by modifying linearly 
the entries of a tensor product of 2 x 2 and 3x3 Hadamard matrices. Such 
constructions result in subfactors having intermediate subfactors, thus of non- 
trivial first relative commutant. Wc analysed computationally our family H{0) 
and we conjecture that it yields subfactors of principal graph A^o- 

In [Ni] the second author introduced a condition for commuting squares, 
called the span condition, which is sufficient to imply isolation of a commut- 
ing square in the class of commuting squares (up to isomorphisms). In par- 
ticular, applying this to commuting squares arising from Hadamard matrices, 
we obtained a sufficient condition for isolation of a Hadamard matrix among 
all Hadamard matrices (up to equivalence). In [TZ] the notion of defect of a 
Hadamard matrix was introduced. Saying that the defect of a matrix is zero is 
equivalent to the span condition. 

It is not settled whether the span condition is also necessary for isolation. 
In [TZ] a possible counter-example is provided: is is shown that no affine family 
stems from the Bjorck-Frobcrg 'cyclic 6 roots' matrix Cg"'' ([Bj],[H]), while its 
defect is non-zero. Therefore, it is asked whether this matrix is isolated among 
all Hadamard matrices. We answer negatively to this question, by showing that 
the family H{6) contains a matrix equivalent to Cg°\ Thus, the question of 
whether isolation is equivalent to the span condition remains open. 

In the last part of the paper we present other new examples of complex 
Hadamard matrices of dimensions 9, 10, 11. These examples were found using 
computers, by doing a numerical search for Hadamard matrices satisfying cer- 
tain symmetry conditions. 

We would like to thank Tcodor Banica, Ingemar Bengtsson, Dietmar Bisch, 
Wes Camp, Romeo Maciuca, Wojciech Tadej and Karol Zyczkowski for fruitful 
discussions and correspondence. Kyle Beauchamp was supported in part by 
NSF under Grant No. DMS 0353640 (REU Grant), and Remus Nicoara was 
supported in part by NSF under Grant No. DMS 0500933. 



2 Self-adjoint Hadamard matrices of order 6 

In this section we classify, up to equivalence, all complex self-adjoint Hadamard 
matrices of order 6. We prove the following theorem, stating that there exists 
a non-afRne one-parameter family of such matrices. 

Theorem 2.1. Let H G Mq{<C) be a self-adjoint Hadamard matrix. Then H 
is equivalent to H{0), for some 6 e [—■ tt, —arcos{ ""^2^ )1 ^ [arcos{ ~^\'^ ),!:], 
where: 
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and the parameters x, y, z, t are given by: 



y = exp{id), z = 



l + 2y-y^ 

yi-l + 2y + y^) 



l + 2y + y^- V2 + 2y + 2y^ + y^ 
l + 2y-y^ 



t : 



1 + 2y + y2 _ v^^l + 2y + 2yi + y^ 



-1 + 2J/ + J/2 

Remark 2.2. In [TZ] it is asked whether the Bjorck-Froberg 'cyclic 6 roots' 
matrix: 
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is isolated among complex Hadamard matrices. It is known that no affine 
Hadamard family stems from C^^"^ ([TZ]). However, this matrix does not satisfy 

the span condition we introduced in [Ni], or equivalently its defect is non-zero, in 

the sense of [TZ]. If cf '^ were isolated, it would follow that the span condition, 
which is sufficient to ensure isolation, is not necessary for isolation. However, 
Theorem 2.1 shows that there exists a continuum of non- equivalent Hadamard 
matrices containing Cq'\ since H{9q) is equivalent to Cq'^ for 9q = 2Arg{d). 
Indeed: 

where P is the permutation matrix: 



P = 



In particular, the question of wether isolation is equivalent to the span condition 
remains open. 

Recall that a complex Hadamard matrix H = (hk.i) G M„(C) is said to be 
dephased or in norm,al form, if hi^k = hk^i = 1 for all k = l,..,n. The next 
lemma shows that, in order to classify all self-adjoint Hadamard matrices, one 
only needs to look at dephased self-adjoint Hadamard matrices. 

Lemma 2.3. Let H e M„(C) be a self-adjoint Hadamard matrix. Then H is 
equivalent to a dephased self-adjoint Hadamard matrix. 
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Proof. Since H = {hk,i)i<k.i<n is hcrmitian, hk.k arc real and thus they belong 
to {—1, 1}. We may assume, by eventually multiplying H by —1, that hi^i = 1. 
Consider the matrix H' = {hk,ihk,ihi,l)i<k,l<n- H' is clearly equivalent to H 
and ^fc 1 = ^'i ; = 1 for all 1 < fc, Z < n. Morover, H' is hermitian since 
^'k I ~ hk,ihk,ihi^i = hi^khi,khi,i = h'l f,. 

□ 

We now recall two easy lemmas involving algebraic manipulations of complex 
numbers. 

Lemma 2.4. Ifx,y,z G C such that \x\ = \y\ = \z\ = 1 andx + y + z = 0, then 
X = ze and y = z^ , where e e + \ - 

Proof. By conjugating x + y + z = Qwe obtain - -\- - + - = Q. Solving and 

X y z 

eliminating x yields — ^ ~ ~ ~^ \- Equivalcntly. (|)^ + | + 1 = 0, which shows 
that y = ze, with e as above. Since l + e + e'^ = Q, x = —y — z = —z{l + e) = 
ze''. □ 

Lemma 2.5. //x, y,z,t £ C such that \x\ = |y| = = \t\ = 1 and x+y+z+t = 
0, then X G {—y, —z, —t}. 

Proof. We have: {x+y){x+z){x+t) = x''{x + y + z + t) + xyz+xyt+xzt+yzt = 
xyz + xyt + xzt + yzt = xyzt{x + y + z + t}=0. □ 

The following lemma is also used in [H], towards the classification of complex 
Hadamard matrices of order 5 . 

Lemma 2.6. Let u,v,s,t be complex numbers on the unit circle. Then: 

{u + v){s + t}{us + vt) e M 

Proof. {u + v){s + t){us + vt) = {us + vt + ui+vs){us + vt) = 2 + {uvst + uvst) + 
{uv + uv) + {st + st) is real, since z + z]s real for every z e C. 

□ 

We now proceed with the proof of Theorem 2.1. Since H is hermitian, its 
diagonal elements belong to {—1,1}. Morover, since for every permutation 

matrix P the matrix H is equivalent to PHP^^, and PHP^^ is still hermitian, 
it is enough to consider the following six possibilities for the diagonal of H: 

Diag{H) E {(1, 1, 1, 1, 1, 1), (1, -1, 1, 1, 1, 1), (1, -1, -1, 1, 1, 1), (1, -1, -1, -1, 1, 1) 
(1,-1,-1,-1,-1,1), (1,-1,-1, -1,-1,-1)}. 

We start by showing that the diagonal of H can not be (1, 1, 1, 1, 1, 1). This is 
the most difficult of the six cases we need to analyse. Indeed, in all the other 
cases the existence of a 1 and a —1 on one of the rows of H will allow us to 
apply Lemma 2.5, thus reducing the number of variables. 
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Lemma 2.7. (a). Let H be a complex 6x6 Hadamard matrix of the form: 
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Then two of x,y,z must be equal. 

(b). Let H e Me{C) be a self-adjoint, dephased, complex Hadamard matrix. 
Then the diagonal of H can not be (1, 1, 1, 1, 1, 1). 

Proof, (a). Assume, by contradiction, x ^ y ^ z ^ x. Denote the last two 
elements on the second and third rows of -H" by u = /i2,5, v = /i2,6j respectively 
s = /i3,5, t = hs^e- Using the orthogonality of the first three rows of H we obtain: 

2 + X + y = —{u + v) 

2 + x + z = -{s + t) 
1 + 2x + yz = — {us + vt) 

Lemma 2.6 implies: 

{2 + x + y){2 + x + z){l + 2x + yz) eR (1) 
The same argument applied to rows 1,2,4 respectively rows 1,3,4, shows: 

{2 + y + z){2 + y + x){l + 2y + zx)GR (2) 

and 

{2 + z + x){2 + z + y){l + 2z + xy)eR (3) 
Expanding the product in (1) and using xx = yy = zz = 1, we obtain: 

{x'^+yz+'l+xy+xz+4x+2y+2z){l+2x+yz) = x'^yz+x'^ + {2xyz+4xyz)+xy+ 
xz+{4xy+xy) + {4xz+xz) + {yz+4yz) + {8x+6x) + {'iy+2y) + {4z+2z) + 13 G M. 

Since {2xyz + 2xyz) + {xy + xy) + {xz + xz) + 4(yz + yz) + Q{x + x) + {2y + 2y) + 
{2z + 2z) + 13 e M, by substracting it from the previous expression it follows: 

x^yz + x^ + 2xyz + xy + xz + 3xy + 3xz — 'dyz + 2x + 2y + 2z & M. 

Thus: 

x'^yz + x'^ +x{x + y + z)- 6yz + {2xyz + 3xy + 3xz + 3yz + 25 + 2y + 2^) e R 
Let 

S = 2xyz + 3xy + 3xz + 3yz + 2x + 2y + 2z 
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Using that {6yz + 6yz) is real yields: 

x'^yz + x'^ + x{x + y + z) + 6yz + S eR 
Similarly, by expanding (2) and reducing real terms we obtain: 

y'^zx + y'^ + y{x + y + z) + 6zx + S eR 

The number S is the same, since (2) is just a circular permutation (.x, y, z) 
{y,z,x) of (1), and the formula for S is invariant to permutations of x,y,z. 
Substracting the two previous expressions and cancelling S, we obtain: 

xyz{x - y) + x'^ - y"^ + {x - y){x + y + z) - 6{x - y)z e R 

Thus: 

{x - y){xyz + X + y + X + y -\- z - Qz) & R 

Hence: 

{x - y){xyz + x + y + zA \ 1 7z) & R 

X y z 

Using that a complex number is real iff it equals its conjugate, we obtain: 



1 1 1 , N /I In/ 1 1 1 1 7, 

{x-y){xyz-\-x-\-y+z-\ 1 1 ~z) = [ )( 1 1 1 \-x+y+z — ) 

X y z X y xyz x y z z 

Multiplying the previous equality by — yields: 

1 1 1 1 1 1 1 7 

-xylxyz + X + y + z -\ 1 1 7z) = 1 1 1 \-x + y + z 

xyz xyz xyz z 

Thus: 

7 , 1 1 1, „ 1111 

- = xyixyz + x + y-\- z^ 1 h-j- "xyz H 1 1 1 \-x + y + z 

z xyz xyz xyz 

Let 

1 111 

T = -7xyz -\ \ \ \ \-x + y + z 

xyz xyz 

111 

it = xyz + x + y + z -\ 1 1 — 

xyz 

We showed: 

- = xyR + T 

z 

By repeating the same argument for the relations (2), (3), and using y ^ z, we 
obtain: 

-=yzR + T 

X 
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where i?, T are the same as in the previous equation, since their formulas are 
symmetric in x, y, z. Substract the last two equations: 

7 7 

= {x- z)yR 

Z X 

Multiplying by -^z^, we obtain: 

7 = xyzR 

This implies = 7. However, by the triangle inequality we have: 

111 

\R\ = \xyz + a; + y + ^;+- + - + -|<l + l + l + l + l + l + l = 7 

X y z 

Since x, y, z were assumed distinct, we can't have equality and thus we have 
reached a contradiction. 

(b). Reasoning by contradiction, assume that H satisfies the hypothesis. De- 
noting H as in part (a) wc know that two of x, y, z are equal. We analyse the 
three possible cases: x = y, y = z, x = z. 
Case I: x = y. Prom (1) we obtain: 

{2 + x + x){2 + x + z){1 + 2x + xz)gR (4) 

If X = y = —1, using the orthogonality of columns 1,2 we obtain /i6,2 = —h^,2- 
Orthogonality of rows 5,6 yields 

2/i6,5 + /i6, 3/15,3 + hQ,ih^,i = 
By Lemma 2.5, we must have 

^16, 5 = ^^6,3^5,3 = —/l6, 4/15,4 

In particular: h^^ih^^z = ft6,4ft6,3- Using this together with the orthogonality of 
columns 3, 4 yields: 

Z-\- /l5,4/l5,3 = 

which together with Lemma 2.4 implies z E {e, e^}, e = exp{2ni/S). It is 
immediate to check that this contradicts equation (3). 

This shows that x ^ —1. Thus, 2 + x + x ^ and after dividing by it in equation 
(4) we obtain: 

{2 + x + z){l + 2x + xz)=4 + 4x + 2x + 2xz + 2xz + x'^z + z eR 

Thus: 2x + x^z + z e K, i.e. ^ + x'^z + ^ = 2x + ^ + z. After multiplying by 
x'^z and simplifying: 

{x^ - l){xz - 1)2 = 

We have x ^ 1, since if a; = y = 1 the orthogonality of columns 1,2 implies 
4 + /i5,2 + /i6,2 — 0, which is impossible. Thus a;2: = 1, so a; = y = ^. It is easy 
to check that in this case relation (3) holds true. 
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Wc substract the sum of the elements of column 4 (which is 0) from the sum of 
the elements of column 3 (which is also 0): 

— /l5,4 — /l6,4 + /l5,3 + ^6,3 = 

Last equation together with Lemma 2.5 yields one of three possibilities: 

I.(i). /i6,4 = — /i5,4 and /le.s = —^5,3- In this case, the orthogonality of 

columns 1 and 3 implies x = —1, which we showed it is not possible. 

I.(ii). ^6,4 = /i5,3 and /i5,4 = /i6,3- Prom the orthogonality of columns 3,4 

we obtain: 



The last equality implies x € {i, —i}. However, the sum of the elements of the 
third column of H is 0: 

2 + 2x + + /l6,3 = 

which contradics the triangle inequality: 2-\/2 = |2 + 2x\ = |/i5,3 + ft6,3| < 2. 

I.(iii). /i6,4 = /J6,3 and ^5,4 = ^5,3- Substracting the inner product of 
columns 4,2 from the inner product of columns 2,3 we obtain: 

^5,2/15,3 — /15, 2/15,3 + ^6,2/i6,3 ~ /i6, 2/16,3 = 

Applying Lemma 2.5, we have three possibilities: 

I.(iii).l. 2/15,3 = /i5,2/i5,3 and /i6,2/''6,3 = /j-e, 2/16,3- Thus /i5,2/i5,3 and 
/i6, 2/16,3 are real. Combining this with the orthogonality of columns 2,3 we 
obtain x real, thus x = ±1, a contradiction. 

I.(iii).2. /15, 2/15,3 = —/i6, 2/16,3- Using this relation together with the orthog- 
onality relation between columns 2,3, we obtain x = —1, contradiction. 

I.(iii).3. /15, 2/15,3 = /i-e, 2/15,3- This equality together with the orthogonahty 
of columns 2,3 implies a; G M, thus x = ±1, contradiction. 
This ends the analysis of the case when x = y. 

Case II: y = z. This case can be treated similarly to Case I. 

Ccise III: X = z. As in the first case, one of the following holds: x = 1, 
X = —1, or X = y = z. However, since the sum of the elements on the third 
row of -ff is we can't have x = 1. Also, if a; = — 1 then equation (2) implies 
y = ±1, which contradicts the orthogonality of rows 2,4. Thus we must have 

X = y = z 

Writing that the sum of the entries of column 2 (which is 0) equals the sum of 
the conjugates of the entries of column 4 (also equal to 0), we obtain: 

/l5,4 + /l6,4 = /l5,2 + /l6,2 

Lemma 2.5 divides now the problem in three cases. 

Case Ill.(i). /i6,4 = — /i5,4 and ft,6,2 ~ ^/i5,2- In this case the orthogonality 
of rows 1,2 forces x = —1, which we showed it is not possible. 
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Case Ill.(ii). ^,5,4 = hf,,2 and ft,6,4 = ^5,2- This case can be treated similarly 
to Case I.(ii), but by looking at columns 2,4 instead of columns 3,4. 

Case III.(iii)./i5 4 = /15 2 and /i6,4 = /i6,2- Again, this case can be treated 
similarly to Case I.(iii), by substracting the inner product of columns 4,3 from 
the inner product of the columns 2,3. 

□ 

Lemma 2.8. Let H e Me{C) be a self-adjoint, dephased, complex Hadamard 
matrix. Then the diagonal of H can not be (1,-1,1,1,1,1). 

Proof. Since the first two columns of H are orthogonal, the sum of the elements 
on the second column of is 0. Two of these elements being —1, 1, the sum 
of the other four equals 0. By applying Lemma 2.4, we may assume, after 
eventually permuting some rows and the corresponding columns of H, that the 
second column of H is (1,-1, a, b, —a, —b). Thus: 
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Using the ortoghonality of columns 3,5 of H, we obtain: 

2y + xa + Z7 = 

Lemma 2.5 implies y — —ax, so z = —xaj. Similarly, considering the or- 
toghonality of columns 4,6 of H, we obtain /3 = —a'f. After making these 
substitutions, the orthogonality of columns 3,4 yields 

l + ab = 

while from the orthogonality of columns 4,5 it follows 

1 - a6 = 

Last two relations are clearly contradictory. □ 

Lemma 2.9. Let H e Mq{C) be a self-adjoint, dephased, complex Hadamard 
matrix. Then the diagonal of H can not be (1, —1,-1, 1, 1, 1). 

Proof. As in the previous lemma, we may assume: 



H = 
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The orthogonality of columns 4, 6, together with Lemma 2.5, implies 

/3 = — a7, z = a'yx 
Since the inner product of columns 3,5 is zero, we have: 

= ax = —y 
The orthogonality of columns 3,4 implies: 

2x = l + db 

and thus ab = x = 1, while the orthogonality of columns 3,6 yields: 

2xa'y = 1 — ab= 

which is a contradiction. 



□ 



Lemma 2.10. Let H G Af6(C) be a self-adjoint, dephased, complex Hadamard 
matrix. If the diagonal of H is (1, — 1, — 1, — 1, 1, 1), then H is equivalent to 
H{6), for some 6 e [-tt, -arcos{ ~^\'^ )\ U [arcos(^i±^),7r], where: 



H{6) = 
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t : 



1 + 2y + - V2v/1 + 2y + 2y3 + y4 
1 + 2y - y2 

1 + 2;/ + t/2 - 72^1 + 2;/ + 2|/3 + j/^ 
-l + 2y + y2 



Proof. By applying lemma 2.5 to the elements on the second column of H, 
which sum up to 0, it follows that the second column of H has to be of one of 
the forms: 
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We may discard the third option, since it is equivalent with the second option 
by permuting rows 5,6 and columns 5,6 of H. Indeed, this operation does not 
change the diagonal of H. By applying Lemma 2.5 to the third column of H, 
which contains a since H is hermitian, we obtain the following possibilites for 

columns 2,3 of H: 
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c 


b 


—a 


b 


— c 


—a 


c 


—a 


—a 


—a 


— c 


-b 




\-b 




\-b 


-a) 


\-b 




\-b 




\-b 


-a) 



We may remove the third arrangement, since it is equivalent to the second by 
permuting rows 5,6 and columns 5,6 of H, and replacing b by —b. We now 
analyse the five cases left. 



Case I. 



A 


1 


1 


1 


1 


1\ 


1 


-1 


a 


—a 


6 


-6 


1 


a 


-1 


—a 


c 


— c 


1 


—a 


—a 


-1 


a 


/3 


1 


b 


c 


a 


1 


7 


V 


~b 


— c 




7 


1/ 



Using that the inner product of the last two columns of if is we obtain: 

-1 + 0.(3 + 2-^ = 

From lemma 2.5 we have: 

7 = 1, a = /3 

Since the sum of the elements of column 5 of is 0, using 7=1 yields: 

3+6+c+a=0 

which implies b = c = a = —1. This contradicts the fact that the sum of the 
elements of column 6is0: 2 — 6 — c + /3 + 7 = 0. 

Case II. 



/I 


1 


1 


1 


1 


1\ 




-1 


a 


—a 


b 


-b 




a 


-1 


c 


—a 


—c 




—a 


c 


-1 


a 






b 


~a 


a 


1 


7 


V 


-b 


—c 


/? 


7 


1/ 



Using the orthogonality of columns 3,4 yields: 

—2c —aQ. — cP = Q 
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Lemma 2.5 implies 

(3 = —1, a ^ —ac 
Similarly, the orthogonality of columns 5, 6 gives: 

7 = —ac 

Using the formulas for a, 7 and the orthogonality of columns 3, 5 we obtain: 

l + ab = 

while the orthogonality of columns 3, 6 yields: 

l-ab = 

which is a contradition. 
Case III. 



H = 



(I 


1 


1 


1 


1 


1\ 


1 


-1 


a 


b 


—a 


-b 


1 


a 


-1 


—a 


c 


—c 


1 


b 


—a 


-1 


a 




1 


—a 


c 


a 


1 


1 


V 


-b 


— c 


P 


7 





The orthogonality of columns 2,4 of H yields: 

26 + aa + 6^ = 
Lemma 2.5 implies aa = bfi = —b. Thus: 

a = —ab,p = —1 

In particular: 

a/? = ab 

However, using orthogonality of columns 1,4 and Lemma 2.5 we have: 
a + /3 + 6- a = 0, thus: {a, /3} = {a, -6} 

which implies: 

a(3 = —ab 

This contradicts a/3 = ab. 
Case IV. 



H = 



(I 


1 


1 


1 


1 


1\ 


1 


-1 


a 


b 


—a 


-6 


1 


a 


-1 


c 


—a 


— c 


1 


b 


c 


-1 


a 




1 


—a 


—a 


a 


1 


7 


V 


-b 


— c 


/3 


7 
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Orthogonality of columns 4,6 yields: 

— 1 — (3 + a'y + p = 0, thus a = 7 
Using this together with the orthogonality of columns 1,5 we obtain: 

2 — a — a = —27 

In particular, 7 has to be real, and since |a + a| < 2 we must have 7 = — 1. 
Thus, a + a = so 

a e {— i, i} 

The sum of the elements of columns 4 and columns 6 is 0. Writing this we 
obtain: 

^ + /3 = 

which shows that /3 e {—i, i}. Using now orthogonality of columns 2,3 of H we 
obtain: 

b = ac 

and writing that the sum of the elements of column 4 is we have: 

ac + c- 0-1 = 
By Lemma 2.5, we have two possibilities: 

ac = 1, (3 = c or c = 1, = ac 

Thus, for each choice of a G {— wc have two possible values of c, /3, which 
uniquely determine the other variables. It is easy to see that all four Butson 
type matrices we obtain satisfy the hypothesis: 



/I 


1 


1 


1 


1 


1\ 


1 


-1 


—i 


1 


i 


-1 


1 


i 


-1 


i 


—i 


—i 


1 


1 




-1 


-1 


i 


1 


—i 


i 


-1 


1 


-1 


V 


-1 


i 


—i 


-1 


1/ 


(I 


1 


1 


1 


1 


1\ 


1 


-1 


i 


1 


— i 


-1 


1 


—i 


-1 


—i 


i 


i 


1 


1 


i 


-1 


-1 


—i 


1 


i 


—i 


-1 


1 


-1 


V 


-1 


—i 


i 


-1 


1/ 


(I 


1 


1 


1 


1 


1 \ 


1 


-1 


—i 


—i 


i 


i 


1 


i 


-1 


1 


—i 


-1 


1 


i 


1 


-1 


-1 


—i 


1 


—i 


i 


-1 


1 


-1 


V 


—i 


-1 


i 


-1 


1/ 



15 



(I 


1 1 


1 


1 


1 \ 


1 


-1 i 


i 


— i 


—i 


1 


-i -1 


1 




-1 


1 


-I 1 


-1 


-1 


i 


1 


i —i 


-1 


1 


-1 


V 


i -1 


— i 


-1 





However, we will see that these matrices are in fact equivalent with a certain 
matrix from the one-parameter family we find in the next case. Thus, it is not 
necessary to include them in the classification. 
Case V. 



H = 



(I 


1 


1 


1 


1 


1 \ 


1 


-1 


a 


b 


—a 


-b 


1 


a 


-1 


c 


— c 


—a 


1 


b 


c 


-1 


a 




1 


—a 


— c 


a 


1 


7 


V 


-b 


—a 


/3 


7 





Wc show that in this case there exists a one-parameter family of solutions. To 
obtain the answer in the form given in the statement of the lemma, let us change 
variables: a = x,b= —y, c = t,'y = z. Thus: 



H = 



(I 


1 


1 


1 


1 


1 \ 


1 


-1 


X 


-y 


— X 


y 


1 


X 


-1 


t 


-t 


— X 


1 


-y 


i 


-1 


a 




1 


—X 


-i 


a 


1 


z 


V 


y 


— X 




z 





Since columns 4,5 are orthogonal, we have yz + (3z = 0, thus: 

(5 = —xyz 

Similarly, the orthogonality of columns 4,6 yields: 

a = xzt 

Using the formula for /3 in the orthogonality of columns 3,6 we obtain: 

1 + xy — xyzt — tz = 
Equivalently: 1 + xy = {xy + l)zi . Using l + xy = l + ^^ = we obtain: 

{xy+l){—-z^={) 
xy 

Wc will assume that xy ^ ±1. Wc treat the case xy — ±1 at the end of the 
proof. Simplifying by {xy + 1) it follows — zi = Q, thus: 

z = xyt 
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Using the orthogonality of columns 1,5 we obtain: 

2-x-t+a+z=0 

and substituting a, z wc have: 2 — a; — t + x{xyt}t + xyt = 0. Equivalently, 
2 — X + y = (1 — xy)t. Since xy 1 we obtain: 

2-x + y 
1-xy 

which implies: 

^ ^ ^~ x + y ^ 2xy + x-y 

Since ti= = 1, it follows: 

2xy + x - y 2-x + y 
xy — 1 i — xy 

Equivalently: 

(2/2 - 22/ - l)x2 + 2(2/2 + 2y + l)a; -{y^ + 2y-l)=0 

Since _ 2y — 1 = does not have solutions of absolute value 1, we must have 
j/2 — 2j/ — 1 ^ 0. Solving the above equation for x wc obtain: 



1 



_ -(2/2 + 2y + l)- V2v/?/ + 2y^ + 2y + l 
y'-2y-l 



_ -(2/2 + 2y+l) + 72^2/^ + 22/3 + 22/ + 1 
2/2- 22/- 1 

where the square root denotes the principal value of the complex power function 

z ^ Z2 . We need to check if these solutions have absolute value 1 when \y\ = 1. 
Consider first the case \xi\ = 1. Denote S = 2{y^ + 2y^ + 2y+ 1). Since \y\ = 1, 
we have: 

5 = 2(l + 2l+2i + l) = l5 

yl y6 y y4 

Thus: 

_ -{y^ + 2y + l)-V~d -(y' + 2y + l) 

1 = XiXi = 



2/2 — 22/ — 1 2/2 — 22/ — 1 

Depending on y, there are two possibilities: 

J -J^ = \^ J = 

V y y V y y 

In the first case, we obtain: 

-(2/2 + 22/+l)-x/^ -{f + 2y+l) + ^^/6 _^ 
2/2 -22/- 1 ■ y2_2y_i 
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which, after substituting 2/ = 1 1 becomes: 

-{y^' + 2y + l)-^^6 -(j/^ + 2t/ + 1) + ^ ^ ^ 
y"^ — 2y — 1 1 — 2y — 

Thus: 

+ 2y + 1) - + 2j/ + 1) + %/^) = (y^ _ 2y - 1)(1 - 2?/ - y^) 

Equivalently: 

+ 2y + 1)2 - 5 = {ij^ - 2y - 1)(1 - 2y - y^) 
It is immediate to check that this identity holds true for every complex number 

y- 

We now show that the case ^ ^5 = leads to a contradiction. We may 

assume S 0, since wc may consider (5 = as part of the first case. By doing a 
similar computation, from |a;i| = 1 we obtain: 

(-(y2 + 2y + 1) - Vd)i-{y^ + 2y + l)-V6) = {y^ -2y- 1)(1 - 2y - y^) 
However, since we showed that 

(-(y2 + 2y + 1) - \/5)(-(y' + 2y + l) + y/6) = {y^ -2y- 1)(1 - 2y - y^) 
for every y, this yields: 

i-iy^+2y+l)-VS){-{y^+2y+l)-V6) = {-{y''+2y+l)-VS){-{y^+2y+l)+VS) 
which implies 

-{y^ +2y+l)-V6 = -{y^ + 2y+l) + V5, thus 6 = 

Cancelation of (— (y^ + 2y + 1) — VS) was possible, since Xi ^ 0. We thus 
obtained a contradiction with the assumption 6^0. 
We have shown that: 



1 if and only if W -^6 = —-^\/S 



y" 



We now need to find for what values of y this holds true. 

Let y = e*^, where Q G (— tt, tt] is the principal value of the argument of y. We 

have: 

5 = 2(y4 + 2y3 + 2y + l) = 2y2((y2 +y-2) +2(y + y-i)) = ^^'\cos{2d) ^2cos(Q)) 
Denote r = cos{29) + 2cos{9). The equality \J -^5 = —^\R> becomes: 
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Denote (j) = Arg{y'^) G (—• tt, tt]. The previous equality becomes: 



If r > 0, we have Vre^ = ^e^'^l'^ and Vre-^ = yfe"*'^/^, thus the equaUty 
we want can not hold. If r < 0, after simplifying by \/--r the equality becomes: 

which clearly holds true. 

We thus need to find the values of 6 such that 

r = 005(26*) + 2005(6*) < 
By denoting p = cos{0), the inequality becomes 

2p2 + 2p - 1 < 

which holds true for p € [ ~^~^ , ~^+^ ]. Since p € [-1, 1], we obtain cos{9) € 
[-l,^i±^], hence 

6 € [— TT, —arcos{ — — ^— — )] U [arcos{ — — ^— — ),7r] 

We have thus obtained a one-parameter family of matrices, which can be easily 
checked to be Hadamard: 

1\ 



H{0) = 



/I 


1 


1 


1 


1 


1 


-1 


1 

X 


-y 


_i 

X 


1 


X 


-1 


t 


-t 


1 


_1 

y 


1 
t 


-I 


1 

y 


1 


— X 


1 

_i 

X 


y 


1 


V 


I 

V 


-t 


z 



1 / 



where: 



y = exp{i9), z 



i + '^y-y^ 



X = 



l + 2y + y^- V2 V^l + 2y + 2y^ + y^ 
1 + 2?/ - 2/2 



t 



1 + 2/y + tj- - y/2 1 + 2/y + 2(y' + 



-l + 2y + y^ 

A similar analysis for \x2\ = 1 leads to another one-parameter family of solu- 
tions, for the same interval of values of 6: 



H'{e) 



/I 


1 


1 


1 


1 


1\ 




-1 


1 

X 


-y 


_ JL 

X 


2/ 




X 


-1 


t 


-t 


—a; 




_ 1 


1 


-1 


1 


1 




V 


t 




y 


t 


V 


— X 

1 
y 


1 

_i 

X 


y 

-t 


1 

z 


1 

1 / 
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where: 



t = 



y = exp{ie), z = /."^.^^ f 2\ 
y{-l + 2y + y^) 

l + 2y + y^ + V2^1 + 2y + 2y^ + y^ 
l + 2y-y2 

1 + 2y + y2 + 72^1 + 2y + 22/3 + y4 



-l + 2y + y^ 

However, it is easy to check that H'{0) is equivalent to H{9): 

PiDiH{e)D2P2 = H'{e) 
where is the unitary diagonal matrix: 




/I 








Vo 



-1 
















1+2 ^1+2 y+2 i 












-y 














1+2 y+y^-V2 sJ\+2y+2 y^+y* 





and D2 is the unitary diagonal matrix: 



/-I 









1+2 y+y'^-V2^1+2 y+2 y^+y* 
-l-2y+y^ 















-1-2 y -y^ + V2^ 1+2 y +2 y'i+y* 
-l-2y+y^ 











yj 










-yJ 



and Pi, P2 are permutation matrices: 



Pi 



fo 


1 











o\ 




fo 








1 





o\ 

















1 




1 





























1 






















1 


1 

















,P2 = 














1 














1 
















1 











vo 





1 












^0 


1 











0/ 



We should also mention here that the four matrices Hi,H2,H3,H4 found in 
case 4, which are easy to check to be equivalent, are in fact equivalent to the 
matrix H{'k/2) which is part of the one-parameter family. Indeed, we have: 



H{'k/2) = P1DH3P2 
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where: 



Pi = 









(I 











\ 
















i 


n 
u 


n 
u 



u 



u 












D = 


- 


-1 



























1 




























- 


-?' 



















\{) 











-V 










(I 





o\ 




















A 





1 















D 


1 














10 






1 























10 


2 = 








1 

















1 












D 





1 





^0 





10 0^ 






^0 


1 











V 



To end the proof, we still need to show that xy = ±1 leads to no solutions. We 
first consider the case xy = —1, thus y = —x and (3 = —xyz = z. Since the sum 
of rows 4,5 is 0, we have: 

2-l+a+a+l+z+z=0 

thus: 

1 + Re{a) = -z 

which implies that z is real, and since \z\ = 1 and |i?e(a)| < 1 we must have 
z = 1, Re{a) = 0. However, since the sum of the elements of row 5 is 0: 

3 = X + 1 — a 

and the triangle inequality implies a = — 1, contradicting Re{a) = 0. 
Consider now the case xy = 1, so y = x and /? = —xyz = —z. Summing up 
columns 5,6 we obtain: 

A^x + t — a — z 

and the triangle inequality shows x = t = 1, a = z = —1. However this implies 
that the sum of elements of column 5 is 1 — 1 — 1 — 1 + 1 — 1 = — 2^0, 
contradiction. 

□ 

Lemma 2.11. Let H £ Me{C) be a self- adjoint, dephased, complex Hadamard 
matrix. Then the diagonal of H can not be (1, — 1, — 1, — 1, — 1, 1). 

Proof. Reasoning as in the previous lemmas, wc may assume: 

/I 1 1 1 1 1_\ 

1 —1 a b —a —b 

^ ^ 1 a -1 X y z_ 

lb X —1 a P 

1 —a y a — 1 7 

\1 -b z P 7 1 / 
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Using the orthogonality of columns 3,5 we have: 

—2y + xa + z'y = Q 

Lemma 2.5 implies 

y = ax, z = 'yax 
Using this and the orthogonality of columns 4,6 we obtain: 

a7 = 

which is not possible since |a| = |7| = 1. 

□ 

Lemma 2.12. Let H G Mq{C) be a self-adjoint, dephased, complex Hadamard 
matrix. Then the diagonal of H can not be (1, — 1, — 1, — 1, — 1, — 1). 

Proof. We may assume 



H = 



(I 


1 


1 


1 


1 


1 \ 


1 


-1 


a 


b 


—a 


-b 


1 


a 


-1 


X 


y 


z 


1 


b 


X 


-1 


a 




1 


—a 


y 


a 


-1 


1 


V 


-6 


z 


P 


7 


-V 



Using the fact that columns 3,5 are orthogonal, as in the previous lemma, we 
obtain: 

y = ax, z = 'yax 

Now the orthogonality of columns 4,6 together with lemma 2.5 yields: 

P = a'y = xz 

Using the expression for /? and the fact that columns 3, 6 are orthogonal we 
obtain: 

1 — ab— z + xfi + yy — z = 
Since xp = z and yy = z, it follows 

1 - aft = 

Thus b = a. But in this case the inner product of columns 3,4 is 

2 - 2x + ya + zl3 = 2 - 2x + 2x = 2 
which contradicts the fact that they are orthogonal. □ 
This ends the proof of Theorem 2.1. 

Remark 2.13. The Butson type matrices Hi, H2, H3, H4,, which are equivalent 
to H{-k/2), are also equivalent to the matrix Dq from the catalogue [TZ]. This 
shows that besides the affine family through Dq, exhibited in [TZ], there also 
exists a one-parameter non- affine family containing Dq. 
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3 More Complex Hadamard Matrices 



We present in this section some new examples of complex Hadamard matrices 
of small orders. These examples were found using Mathematica, by searching 
for the local minimum of a function f{H) encoding the conditions a matrix H 
needs to satisfy to be Hadamard. i.e. f{H) = if and only if H is Hadamard. 

Trying to find such general solutions leads to numerous numerical results 
that are difficult to interpret formally. However, it turns out that asking for 
extra symmetries for H, such as H hermitian, often yields some clear algebraic 
results. 

Inspired by our work on the classification of 6 x 6 Hadamard matrices, we 
tried to fix certain diagonal entries of H, making them 1 or —1. Such con- 
ditions seem to be very strong, leading to solutions that are easy to interpret 
algebraically. 

New examples for n=9. 



Ha = 



(I 


1 


1 


1 


1 


1 


1 


1 


1 \ 


1 


-1 


-1 


y 


y-' 




y-' 


y 


y-' 


1 


y 


y3 


-y' 




y 


-1 


-y 


y 


1 


y 


y 


-1 


y-' 


y-' 




-1 


y-' 


1 


-y 


-y' 


y' 


y' 


-1 


y 


y 


y 


1 


y 


y3 


y' 


y 


y-' 


y-' 


y 


y' 


1 


y' 


y 


y 


y' 


y-' 


y-' 


y' 


y 


1 


-y' 


-y 


y 


y 


-1 


y 


y' 


y' 


V 


y' 


y 


-y 


y 


y 


-1 


-y' 


y"" ) 



y = 



_ l-iVl5 



It can be easily checked that Hq satisfies the "span condition" of [Ni], thus 
it is isolated among Hadamard matrices. In particular, it is not equivalent to 

any of the matrices from the 4-parameter family Fg^'' described in [TZ]. Also, 
Hg is not equivalent to any of the known circulant solutions, since its entries 
belong to 

The next matrix was obtained by searching for the local minimum with a 
fixed diagonal: 



BNo 



(I 


1 


1 


1 


1 


1 


1 


1 


1\ 


1 


-1 




e3 


-1 




e« 




e 


1 




-1 




e 


£3 


-1 






1 


e3 




-1 


e 


e8 


e9 


e3 


-1 


1 




e 


-1 


-1 


e3 




£2 




1 




-1 


e 


£3 


-1 


e 




e6 


1 


e 




e9 




e 


-1 


-1 


e3 


1 










-1 




-1 


e 


^1 


-1 


£2 








e3 


e 


-y 



g _ g27ri/10 



The matrix BNg has defect 2, in the sense of [TZ], thus it might be part of 
a family of non-equivalent Hadamard matrices. 
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New exEimples for n=10. 



BNw 



(I 


1 


1 


1 


1 


1 


1 


1 


1 


1 \ 


1 


1 


a-2 


a-i 




a2 


1 


a2 


1 


a 


1 


a" 


1 


a-2 




1 


a 


a2 


1 


a-i 


1 


1 


1 


1 




a2 


a2 


a 




a-i 


1 


a2 


a-2 


1 


1 


a 


1 






a-i 


1 


1 


a-i 


a-2 


1 


1 




a2 




a 


1 


a2 


1 


a-2 


a-i 


a2 


1 


1 




n 


1 




a-2 


1 


a-2 


1 




1 




a 


1 


a 


a-2 


a-2 


1 


a2 


a2 


1 


1 


a-i 


^1 


a 


a-i 


a-i 


a-i 


a 


o 


a 




1 / 



a = 



_ -l+»v/l5 



This matrix was found by numerical search of the local minimum, with the 
constraint that all diagonal entries be 1. It satisfies the span condition and thus 
it is isolated among complex Hadamard matrices. 

New exeimples for n=ll. 



( 1 


1 


1 


1 


1 


1 


y 


y 


a 


b 


1 


y 


b 


b 


a 


1 


a 


b 


y 


y 


1 


b 


a 


y 


b 


1 


b 


y 


b 


y 


1 


-1 


X 


-1 


X 


1 


X 


X 


-1 


-1 


1 


-1 


-1 


X 


-1 


1 


X 


-1 


-1 


-1 


1 1 


-1 


-1 


X 


X 



where x 
are isolated. 



|-i^and 











or 




-X-) 







1 


1 


1 


1 


1 


\ 


-1 


X 


-1 


X 


-1 




X 


X 


-1 


-1 


-1 




-1 


-1 


X 


-1 


X 




X 


-1 


-1 


-1 


X 




-1 


-1 


X 


X 


-1 




X 


X 


X 


X 


X 




a 

X 


b 

X 


y 

X 


y 

X 


b 

X 




b 

X 


y 

X 


b 

X 


y 

X 


a 

X 




y 

X 


b 

X 


b 

X 


a 

X 


y 

X 




y 

X 


y 

X 


a 

X 


b 

X 


b 

X 


) 


b 


a 


y 


b 


y 



1 



- 1 . Both these matrices 
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